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Abstract 

We propose a general scheme of constructing of soliton hierarchies from finite dimensional Stackel 
systems and related separation relations. In particular, we concentrate on the simplest class of 
separation relations, called Benenti class, i.e. certain Stackel systems with quadratic in momenta 
integrals of motion. 

> . 

CN ■ 1 Introduction 

The theory of integrable nonlinear evolution equations has a long history as a part of many branches 
of theoretical physics and applied mathematics. Generally it can be divided in two parts: the theory 
jy-j , of integrable nonlinear ordinary differential equations (ODE's) and the theory of integrable nonlinear 

partial differential equations (PDE's). Within the first class of equations (ODE's) we will consider finite 
dimensional Hamiltonian systems, integrable by the Hamilton-Jacobi method, called Stackel systems, 
while within the second class (PDE's) we will consider (l-fl)-dimensional field systems, having infinite 
hierarchy of commuting symmetries and called further for simplicity soliton systems. The solvability by 
quadratures of some class of finite dimensional systems by the Hamilton-Jacobi method, laid in the 19-th 
century one of the fundaments of analytical mechanics of integrable systems, while the solvability by 
quadratures of some class of infinite dimensional field systems by the Inverse Scattering Method, laid in 
second half of the 20-th century one of the fundaments of the so called soliton theory. 

During the last few decades both theories have been developed very intensively using many common 
modern mathematical tools like Lax representation, r-matrix theory, multi-Hamiltonian theory etc. In 
that time some links between both theories were investigated. It was found (P^-U, see also references 
in |S]) that finite dimensional restrictions, invariant with respect to the action of a given soliton system, 
like stationary flows, restricted flows or constrained flows of Lax representation, are Liouville integrable 
Hamiltonian systems of Stackel type. Moreover, analytical solutions of an appropriate finite dimensional 
systems are closely related to a special class of solutions of related soliton systems, like for example 
so-called finite-gap solutions HI],[7|. 

In the present paper we are interested in passing in the opposite direction - building integrable 
hierarchies of PDE's from Stackel systems 8 . In that sense we would like to initiate a unified approach 
to Stackel ODE's and soliton PDE's. Our claim is the following: both a wide class of Stackel systems and 

t*Partially supported by KBN Research Project 1 P03B 111 27 and by Swedish Research Science Council grant no 
2004-6920. 
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a wide class of soliton systems can be constructed form common fundamental objects known as separation 
relations (or from separation curves). 

The paper is organized as follows. Section 2 is devoted to general description of the concept of 
separation relations. Section 3 explains the main ideas of relating soliton systems with separation curves 
that are quadratic in momenta. The idea is to apply to a set of Killing vector fields a set of invariants 
generated by Euler-Lagrange equations associated with appropriately chosen Lagrangian densities. This 
allows for elimination of some variables in our Killing systems which leads to dispersive soliton hierarchies. 
Section 4 is a brief introduction to what can be called Benenti class of Stackel systems. In section 5 we 
describe the structure of our systems in Viete coordinates. In Section 6 we explain the details of our 
elimination procedure which allows, in a systematic way, to construct soliton hierarchies. It is divided 
into two subsections as the elimination procedure differs in case of " positive" and " negative" (see below) 
separable potentials. Section 7 concludes the article with several examples. 



2 Separation relations 

Let us consider a 2n-dimensional manifold M equipped with a Poisson operator II with some canonical 
(Darboux) coordinates labelled as M 3 u = (//, A), with fi = (/ii, . . . , [x n ) and A = (Ai, . . . , A n ). The 
following definition introduces the basic object of our considerations UJ: 

Definition 1 A set of n relations of the form 

iPiiXufafCtu. ..,a n ) = 0, i = 1, .. .,n, a, e R (1) 

(each involving only one pair Ai, jii of canonical coordinates) are called separation relations provided that 
the dependence of ip on a is essential i.e. that det 

The condition in Q means that we can resolve the equations with respect to a.; obtaining 
ai = Hi(X,fi), i = l,...,n. This defines some new functions Hi(X, /i) that in turn define the follow- 
ing Hamiltonian systems (evolutionary vector fields) on M 

u u =ndH i = X Hi , i = l,...,n. (2) 
If the functions Wi(Xi, a) are solutions of a system of n decuple ODE's obtained from JTJ by substituting 

_ 3ff,(A„q) 
I 1 ' 1 ~ dXt 

(fi\Ai,fii= 7^- ,oi, . . . ,a n I = 0, i =l, ...,n, (3) 

then the function W(X, a) — X)™=i^ / *(^' a ) ^ s an additively separable solution of all the equations J2J 
and simultaneously it is a solution of all Hamilton- Jacobi equations 

a,=^(A, M = «M), i = 1 ,..., n (4) 

related with the Hamiltonians Hi - simply because solving Q to the form a, = Hi(X, n) is a purely 
algebraic operation. Assume now that det f gXi&a ) = ( d\ da ) ^ ^" Then the Hamiltonians Hi 
Poisson-commute since the constructed function VK(A, a) is a generating function for the canonical trans- 
formation (A, fi) — > (6, a) where 

dW(X,a) 

bi = —tt + constt, i = ±, ...,n. (5) 

oai 

Equations JSJ) are implicit solutions of J5| known as the inversion Jacobi problem. Thus, starting from 
a set of n separation relations we can create an n-dimensional separable Liouville system. All systems 
separable in the sense of Hamilton-Jacobi theory can be obtained in this way. 
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In an important case, when the functions tpi in Q do not depend on the index i, the separation 
relations Q can be generated by taking n copies of a curve in A-/i plane: 

ip(X,IJ,,ai,...,a n ) =0, en € R (6) 

called separation curve. 

Restricting our considerations to a subclass of Jy, when all separation relations are afhne in <Zj = 
with coefficients being monomials in A and /z, we obtain 

n 

Y,H k tf k \? h =^(Aj,ft), « = l,...,n, a fe ,/? fc eN (7) 
fe=i 

where ^ are arbitrary smooth functions of two arguments. Equations Q are called generalized Stackel 
separation relations and the related dynamic systems, generated by Hamiltonian functions Hi, are called 
the Stackel separable ones. To recover explicit Stackel form of Hamiltonians it is sufficient to solve the 
linear system Q with respect to Hi. If additionally ^,-(Aj,//j) = ^(Aj,/Uj) then the above separation 
conditions can be represented by n copies of the following separation curve: 

n 

Y^Hu^X^ =V(A,M)- (8) 
fc=i 

The separable systems that were most intensively studied in the last century were one-particle dynamical 
systems on Riemannian manifolds with flat or constant curvature metrics. All these systems can be 
obtained by choosing on = 0, = n — i, i = 1, n with ip quadratic in momenta 

^(A, Ai ) = i/(A) M 2 + 7 (A). (9) 

This case will be considered in the next sections of this article. 

We can now shortly present - by a simple example - the possibility of passing from a separation curve 
to soliton systems JH]- Let us consider the separation curve (jHJl with n = 2, ax = oli = 0, (3% = 1, (3% = Q 
and with ip in the form of 

H 1 X + H 2 = ^A^ + A 4 . (10) 

The related separation conditions J7J) are 

HxXi + H 2 = oAiA*! + X\ , . 

ffiA 2 + ff 2 = |A 2 /i| + A| 

Solving this linear system with respect to Hi and £f 2 one gets the Liouville integrable system (j2J on four 
dimensional phase space, written in separation coordinates (A,/i). The explicit form of Hamiltonians Hi 

1 A lM ? - X 2 (i 2 2 + 2A? - 2A| „ 1 AiA 2 (/x? - t4 + 2A? - 2A|) 
-"l = t: ; ; 1 "2 _ 



2 Ai — A 2 2 A 2 — Ai 

The canonical transformation of the form 

1 , 

gi = Ai + A 2 , ~q 2 = -AiA 2 , 
AiMi A 2 ^i 2 ; — — — / ^ \x 2 



Pi = \ r- + i r") P2 



^^(a7 



Ai — A 2 A 2 — Ai ' \Xi — A 2 A 2 — Ai 

transforms the system to new coordinates (q,p), with 

Hi = ^p? + ip 2 + ql + -^qiql, H 2 = ^q 2 pip 2 - \qiPl + JqQ2 + j1iQ2- 
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The function Hi(q,p) turns out to be the Hamiltonian function of the integrable case of the Henon-Hcilcs 
system, while H2(q,p) is the additional involutive first integral of this system. Let us now denote the 
evolution parameters t\ and ti of the system by x and t, respectively. Then we obtain 

Qx,x = p — =Pi, 92,x = -a — =P2, 

<9£T 2 1 <9# 2 1 

qi,t = — = 7T<?2P2, Q2,t = — = o<72Pl - q\P2, 
dpi 2 dp 2 2 

from which eliminating p\ and p 2 we obtain a system of first order PDE's for qi(x, t) and q2(x, t) 

Ui,t = ^92<72« = tC^Dx, 32,* = 2?29i,x - gi^.x- (12) 
Finally, we can eliminate 92 through 

dHi . 2 1 2 
ogi 2 

which yields q\ = —Qqf — 2qi, X x and then generate a higher order (in a;— derivatives) PDE. The first 
equation in (|12|l turns then into the famous KdV soliton system 

qi.t + ^qi,xxx + 3qiqi,x = 0, (13) 

while the second equation in (|12|l turns into a differential consequence of the first one. Obviously, just 
from the presented construction, there is no guarantee that equation lfO|l is integrable. We can only say 
that qi(x,t) calculated from the corresponding inverse Jacobi problem is a nontrivial particular solution 
(one-gap solution) for the field system Ijl3|l . To prove the integrability of (|13|l one has to construct some 
related infinite hierarchy of symmetries using some more regular procedure. 



3 From separation curves to constrained dispersionless systems 

In this paper we will concentrate on a special but important class of separation curves with the function 
ip(X, n) in JSJ being quadratic in momenta /x, (more precisely, of the form and with multipliers of 
Hamiltonian functions being monomials with respect to A 

H 1 \^+... + H n A' 3 " = i A' V + A fc , (14) 

where (3\ > ... > /3 n _i > j3 n = 0, Pi £ N, m, k £ Z and n £ N. Separable systems from this class 
describe one-particle dynamics on Riemannian manifolds and belong to classical Stackel systems. Each 
class of these systems is labelled by a decreasing sequence (/3i, (3 n ) while members of a given class are 
numbered by pairs (m, k) £ Z 2 . Taking n copies of the curve (fT4*j) with variables (A, n) labelled within 
each copy as (A,,/Zj), we obtain a system of n separation relations in the form of n equations linear in 
the coefficients Hi. Solving it we obtain n functions Hr m = Hf 71 (A, fx) of the form 

H (m,k) = - fl T KrG (m) fx + y(k) t r = l j ... jnj m 7 k£Z (15) 

where we denote A = (Ai, . . . , A„) T and fi — (jtii, . . . ,/i„) T . The functions l|15|l can be interpreted as 
n Hamiltonians on the phase space T*Q cotangent to a Riemannian manifold Q equipped with the 
contravariant metric tensor G^ m \ These Hamiltonians are in involution with respect to the canonical 
Poisson bracket on T*Q. Moreover, they are separable in the sense of Hamilton- Jacobi theory since they 
by the very definition satisfy Stackel relations 114(1. The objects K r in (|15|l can be interpreted as (1, 1)- 
type Killing tensors on Q. The scalar functions Vr k ^ are separable potentials. Further, all the metric 
tensors G^" 1 ' and all the Killing tensors K r are diagonal in A- variables so that: 

tf r = diag(^,...,0 (16) 
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where v\, are eigenvalues of K r . We will constantly assume that these eigenvalues are single. 
The set (|15f) of n Hamiltonian functions leads to n Hamiltonian systems on T* Q of the form 

K = ^H~' ^ = BX~> r = l,...,n. (17) 

Let us now call the variable t\ as x; t\ = x. Since all the Hamiltonians Hr k,m ^ (for fixed k and to) 
commute, the equations (|17|l have a common set of solutions depending on all the evolution parameters 
U 

Aj = \i(t\ = X, <2, ■ • ■ , t n ), Hi — Hi(ti = X, t2, ■ ■ ■ , t n ). 

We have, due to (fTTft . that 

f)Tr(k,m) 

X x = A tl = — i = G< m V so that h = g^\ x , 

o/i 

where the inverse of (i.e. respective covariant metric tensors) is denoted as gt m \ Observe that the 
above relation does not depend on k. Using this to eliminate /z from the first part of (|17fl we obtain 

oo-(fe,-m) 

on 



or, according to the above 



X u =K r X x = Z?(X,X x ), r = l,...,n. (18) 



This is a set of n autonomous systems of n coupled first order PDE's of evolutionary type, with the right 
hand sides depending linearly on the derivatives X x More precisely, it is a set of n integrable disper- 
sionless equations, belonging to the class of so-called weakly nonlinear semi-Hamiltonian systems |12 | .|13 | . 
where the variables Xi are the Riemann invariants for (|18|) . We will call them Killing dispersionless system 
as they are constructed directly from Killing tensors. 

We will interpret the right hand sides of (|18|l as vector fields on an infinite dimensional manifold M. 
the points of which are vector functions of x of the form u = (Ax {x), . . . , A„(a;)), where we assume that the 
functions Xi(x) are either periodic in x or they vanish together with all their derivatives when x — » ±00. 
A vector field X is at a point u £ M given by an n-tuple of the form X(u) = (/i[A], . . . , / n [A]) where 
fi[X] — fi(Xi, Xi. x , . . . , A2, X2, x , ■ ■ ■ , X n , X n . x , . . .) are differential functions of A. Similarly, a covector field 
a on M. is in a point u = (Ai(ie), . . . , A„(a;)) given by a(u) = (gx [A], ... , g n [A]). The dual map between 
T U M and T*M is given by 



J X 



Here and below the integration is performed over one period (in case of periodic boundary conditions) 
or over R in case of functions vanishing at ±00. All functions and expressions are always assumed to be 
integrable. For any two given vector fields X and Y on Ai their commutator is defined in a usual way as 
[X, Y] — X' [Y] — Y' [X] where X' [Y] denotes the directional derivative of X in the direction of Y. 
As was shown in the vector fields Zf pairwise commute: 

[Z?,Z?]=0 i,j = l,...,n, 

thus, (|18|l is a set of n commuting evolutionary dynamic systems (vector fields) on M.. We will need 
the superscript n to indicate the number of components (dimension) of these systems. Below we will 
introduce invariants on i|18|) that eventually turn these systems into hierarchies of soliton systems with 
lower number of fields. This is the main idea of this paper. 

We begin by defining the following differential functions (currents, 'Lagrangians'): 

4"' m ' fc) = f \x T x g^K r X x - V} k \ r = 1, . . . ,n. (19) 
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In our further considerations we will especially need the first current £^ n ' m ' k ^ so we w [\i denote it simply 

by£ (n,m,fc) ; 

£(n,m,k) dof ^,(n,m,fc) 

This current is a Legendre transform of H[ n ' m ' k ^ (this is not true for _£^™> m >' £ ) with r > 1). These 
differential functions yield the following functionals on M.: 

li n ' m > k) {u) d = f £l n ' m > k) [\]dx, 

J X 

where, as usual, u = (Ai(x), . . . , X n (x)), We have, of course, 



dt s \ 5\ 



where £ = (Ei, ...,E n ) = yj^, ■ ■ ■ , -gj-J * s ^ e Euler-Lagrange operator on A4. 
Lemma 2 In the notation as above, 

- k ^—=Q, r = l,...,n. 

Proof. It suffices to prove that 



n n 

Ei (c^ m ' k ^ \ lAr =J2 E i (4"' m ' fc) ) (20) 



; = 1 



since integrating of (|2Tj|) yields, 



dt r dt\ 

while 



d An,m,k) 



III 7 n } n p J 

— = j SUE, (4»->) A,., * = I - (£[•-») * = 0. 



due to the appropriately chosen boundary conditions. The proof of (|20|l can be found in Appendix A. ■ 

Corollary 3 Lemma\^ due to theorem of I15f (see also )lftf ) implies that the 2n- dimensional set £ C M 
defined as 

£ = {u: Ei (£(™< m ' fe )) = for all i = 1, . . . , n\ 

is Z™ -invariant for all r — 1, . . . , n. 

Thus, if Uq G £ then the integral (Frobenius) n-dimensional submanifold S uo of M. spanned by the 
commuting vector fields and containing m is a subset off. This means that the solution X(x,t r ) of the 
r-th Killing system in l|18|) that starts at a point uq £ £, i.e. initially satisfying the set of Euler-Lagrange 
equations 

Ek U( n < m > k A =0, i = l,...,n (21) 

remains in £ . i.e. always satisfy 121|) . This further means that we can use the set of equations l|21|l to 
eliminate some of the variables in (|18fl . Such an operation does not alter IjlSjl . but reparametrizes 
it, leading to fewer equations of higher order, and the dispersion will occur. As we will see below, this 
operation of elimination of variables from (|18fl through the use of (|21|) will lead both to known and new 
soliton hierarchies in (1 + 1) dimensions. 
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4 Benenti class of Stackel systems 

In the rest of this paper we consider the simplest class of separation curves 1(14(1 in the form 

ffiA"- 1 + H 2 X n - 2 + -..+H n = ^ A "V + ^ (22) 

(A, /i G R for a moment), where n € N while m, k £ Z. This object contains a complete information 
about the so-called Benenti systems ^7J-Qni- Hamiltonian functions calculated from the related system 
of separation relations take the form 1(15(1 f2U] . Due to a special form of 122(1 it turns out that the metric 
tensors G^ are now 

G (m) = £ m G (0) ) with G (0) = diag -L\ : m eZ, 

\Ai A n J 

where Aj = Yl (\ ~ A?) ancl where L =diag(Ai, . . . , A n ) is a (1, l)-tensor on Q (it is a conformal Killing 

tensor with respect to all the metrics G^). Moreover, Killing tensors K r can now be obtained by the 
following recursion relation: 

K r+1 = LK r + q r I, Ki=I, K n+1 =0,r=l,...,n, (23) 

so that indeed they are diagonal (in A-coordinates) in accordance with 1(16(1 : K r =diag(i^, . . . , V™). The 
functions q r = q r {\) are coefficients of the characteristic polynomial of the tensor L i.e. they are defined 

by 

n 

det(AJ-L) = ^2 qi X n -\ (24) 

so that qo — l,q± — — Y^i=i ■ ■ • > In — (~ l)"Ili=i^ (?* are Viete polynomials in the variables A). 
Moreover, the potentials in Hamiltonians 1(15(1 can now be obtained from the following recursion 
relation |2()| : 

7W = y r ( |- 1) - gr yf- 1) , fcez (25) 

(with the convention that Vr ^ = for r < 1 or r > n ) with the initial condition: 

K (0) = $r,n, r=l,...,n. (26) 
This recursion can be reversed. The inverse recursion is given by 

V W = ~ —V 7 [ k+1 \ /,- , Z r ! n. (27) 



The first potentials arc rather trivial: 



K (fc) = Sr,n-k for fc = 0, 1, n - 1, K (n) = ~9r, = (28) 

On 

but for r < — 1 or for r > n the potentials become complicated polynomial (for r > n) or rational (for 
r < —1) functions of q. 
From 1)25(1 we get 

V„ ( * ) =-^ 1 ( *- 1) J fcez (29) 

and 

K (fe) = -g,^ fc - 1) -g r+1 ^ fc - 2) -...- g „y 1 (fc -" +r - 1) , fcez (30) 

while iteration of (|2"7jl leads to 

= -- U-iV}^ + ... + + v,[ k+r A 

= q r -iv} k) + ... + qi v 1 {k+r - 2) + v} k+r ' 1 \ kez. (31) 
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5 Killing systems and related invariants for the Benenti class in 
Viete coordinates 



The functions <2v(A) defined in H24|) can serve as a new set of variables on Q (we will call them Viete co- 
ordinates), ft turns out that these coordinates (that also reparametrize the infinite-dimensional manifold 
M so that M 3 u — (qi(x), . . . ,q n (x)) now) are much more convenient for our further purposes. The 
above considerations, in particular Lemma[21and the corollary that follows, remain true independently of 
coordinate system since the Euler-Lagrange equations are invariant with respect to point transformations. 
In this section we sort out the structure of i|18|) and (|21[1 for the Benenti class in Viete coordinates as 
well as we prove many other important relations. 

The tensors L, L~ x , G^and g^ — (G^ ^) have in Viete coordinates l|24|) the form: 



L 



G (o) 



so that L\ = V> 



( -qi i 

-<Z2 o 

V ~<ln 



/ 

••• 

1 
V 1 91 

and gf = V[ 



1 

o / 

1 \ 



/ ••• 
10 



L~ = 



,(o) 



Qn-l j 



V o 

/ vl 2n - 2) 
V 



i \ 
) 



- ILL 



(32) 



-91 1 \ 
•• 



-qi 
l 



(33) 



J 



(2n-i-j) 



Viete coordinates the form q tr = K r (q)q x or, explicitly 



Moreover, for the Benenti class, the system (|18(l attains in 



dt 



= (qj+r-l )x + ilk (qj +r -k-l) x - Qj+r-k-1 (%) x ) = (K [q]Y T, j = 1, . 



(34) 



fe=l 



where q a — as soon as a > n and (Z™ [q]Y denotes the j-th component of the vector field Z™ [q\. One 
proves l|34|) by a direct calculation, using (|32fl and l|23|l . Observe, that the following symmetry relation 
takes place: {Z? [q]) ] = {Z™ [q])\ i,j = 1, . 



,,n. 



We can, in accordance with the above, also define the infinite Killing hierarchy for the Benenti class 



-jj-qj = (qj+r-i)x + Y (qj+r-k-i) x - qj+ r -k-i (%) x ) = {z? [qJY r,j = 1,. ..,00 



(35) 



fc=i 



that is formally given by the same expression as Ij34(l but where we now do not impose the restriction 
q a = for a > n. By comparing (|34|l and (|35|l one sees directly that for the r-th Killing vector field 
Z™ [q] from (|34|l its first n + 1 — r components coincide with the corresponding components of the infinite 
vector field Z™ [q]: 



(Z? [q]Y = {Z™ [q)Y forj+r-l< 



(36) 



Lemma 4 The infinite- component vector fields Z^ 2 [q] in mutually commute: 

[Z°° [q] , Zf [q]] — for all i,j — 1, ... 00 

Proof. This can be proved by using [Zf [q] , Zf [q]] = for alH, j = 1, . . . n and Indeed, from © 
and the relation 

(Z?°[q}) j = (Z°°) j [qi,-,qi +j -i] 
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one finds that 

( M , Zf [q]) ) ' = ( [Zf n ~ 1] [q] , Zf n ~ 1] [g]] ) ' , i, j, I = 1, .., n 
for arbitrary n 6 N. ■ 

Let us point out that the infinite Killing hierarchy Ij35ll is exactly the so-called universal hierarchy 
considered recently in |21).|22| from the point of view of Lax representation. 

Lemma 5 In Viete coordinates the following relations hold: 

1. 

dV, W dv! k+a) 

— 1— = for i = 1, . . . , n - a, k e Z (37) 

c% aq l+a 

2. 

(L k )]=vt n+k - 3 \ fceZ (38) 

3. 

g {m) = y{2n-m-i-j) ^ mgZ (3g) 

Proof. For relation 13 7|) the proof is inductive with the help of formula l|30|l . For relation l|38[) the proof 
is by induction with respect to k. By L l - — By the induction assumption and due to the 

recursion (|25|l 

n 

(L k+1 )] =J2 L r (L% = -qiV} n+k - j) + V^ k - j) = Vl {n+k - J+1) 

r=l 

which concludes the inductive step up. Similarly, due to the recursion 127|) . 

n 



r=l 



which concludes the inductive step down. Finally, for relation l|39|l . according to (|33|) . we have g, 
y( 2n 1 i) _ gy induction 



(o) 



4 m+1) = Y j Vl 2n - m - i ~ k) {L- 1 )] 



fe=i 

Thus, due to (|3*2|) we have for j < n 



fm+l) (m) T r (2ri— TTt— i— j— X) 

9ij -«j+i = V? 



while for j = ?i we have 
(m+i) 



1 / Tr (n — rn— i) . T7 -(2n— m— 2— i) , T r (2n — rn— 1 — i)\ 

— lg n -iVi +... + ' + I 

qn ^ ' 

_y(n— m— i) _ y(r»— m— 

yn 

which follows from (|29|) and l|31(l . This concludes the inductive step up. The induction down (for m < 0) 
is proved in a similar way. ■ 

The next theorem describes symmetry properties of functions Ijl9l) . Observe that due to l|28|l the 
functions (|19fl are in the Benenti case geodesic (without the potential part) for k = 0, . . . , n — 1. 

Theorem 6 For the Lagrangian densities 

-i n 

Ak) 



the following relations hold: 
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1 . for a = 1, . . . , n — 1 



Ei (£ n < m < k ) = E^ a (£ n > m+a < k - a ) , i = a + l,...,n, (40) 
that can also be written as 

E . (£n,m,fc) = E . +a (£n,m-a,fc+ Q ) ^ • _ ^ _ >n _ a . ( 41 ) 

2. 

(£«.°.2n+^ = ( jC „+l,0,2n+a+2^ (7 = 1, 71 — 1, I = a + 1, n. (42) 

3. 

£ ; (£"<"- CT >°) = (£»+i.«+i-».o) j ,7=1, ...,„_!, ? = 1,...,ct (43) 



(£"."-■- ")|^ +1 = £7, (£«+i.«+i-,-«-i) , / = 1, n -a, j= 1, n (44) 

The proof of this theorem can be found in Appendix B. As it will be shown in the next section, 
Theorem [5] guarantee that the form of invariants survives the passage from n-component to (n + 1)- 
component Killing system and hence it will be crucial for the construction of soliton hierarchies. The 
index a will be related with the number of components of the obtained soliton systems. 

6 Elimination procedure 

Using the results of the previous section we will now construct in a systematic way soliton hierarchies 
related to Benenti class of separation relations. These hierarchies will be generated by a procedure of 
elimination of variables in the set of dispersionless systems (|34() with the help of Euler-Lagrange equations 
(|21|l (with suitable chosen parameters n, m and k determining the metric tensor g( m ) and the separable 
potential v/ fe ^). Actually, we present two separate elimination procedures, one for positive potentials 
(i.e. with k > n) and one for negative potentials (i.e. those with k < 0), leading to different soliton 
hierarchies. 

As we pointed out, the set £ C M. of solutions of Euler-Lagrange equations (|21|l is invariant with 
respect to all the vector fields Z r of Killing systems l|18(l . The same must be true even in Viete coordinates: 
(|21|l written in Viete coordinates is also invariant with respect to Ij34(l . This means that we can use the 
set of equations l|21|l to eliminate some of the variables qi in 1)34(1 . since along the solutions of systems 
from (|34(l they are all the time mutually related by the relations (|21|l . 

6.1 Elimination for positive potentials 

First, let us concentrate on the case with positive (polynomial) separable potentials. Our aim is to 
produce s (s S N) commuting cr-component [a £ N) vector fields (evolutionary systems) from l|34|l and 
l)2ip. In order to do this we choose nasn = s + a — 1 and consider the set of systems (I34|l with this 
chosen n: 

q tr =Z? + °- 1 [q li ...,q s+a _ 1 ], r = 1, . . . , n = s + a - 1. (45) 

Notice that r can reach r = s + a — 1 but only up to r = s the first a components of (|45|) are complete 
in the sense of the infinite Killing hierarchy (|35(l . i.e. coincide with the corresponding components of this 
hierarchy (see (jSHJ). We will now use 121JI generated by £ n >°. 2,l +°' j n order to perform the elimination. 
The structure of these equations is described in the lemma below. 

Lemma 7 The last n — a invariant equations in f21\) for £™-A 2 ™+ CT j with n = s + a — 1 (so that m = 
and & = 2n + (7 = 2s + 3<7— 2) have the form 
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(n.a) r, (n,<r) r 1 n 

K+i = - 2 <in + Paii'm, -,q n -i\ = o 

: (46) 
w' n ' = -2q a+1 + tp^n '[qi,...,q a \=U 

where we denoted, to shorten the notation, Ei (£ n ,°2 n +°' s ) a s wj n ' a ^ . 
Proof. From the recursion 125(1 it follows that 

v ( n+j) = v (n +j ) (qu . = Q) ^ _ x (4?) 

so that, again by (pT7|) and (|SU|) 

Qy(2n+a) Qy( n +3) 

— ^ = — i = 2q j + f j (q 1 ,...,q j _ 1 ) j = 2,...,n, a = 1, ...,n- 1 (48) 

oq n +i-j+a Oqi 

(for j ' = 1 we would have /i = 0) where the first equality follows by inserting i = l, k = n + j and 
a = n — j + er in (|37|l and the second one from the fact that according to ((28(1 and l(3U(l V x +,? ' = 
—qj+i + 2gi gj + tfj(qi, <Zj-i) for j = 2, .. . ,n. On the other hand, for the geodesic Lagrangian density 

n 

*- — 2 / y K i qi,xqj,x 

from (E2J, as v[ k) ^ for fc > n- 1 and dV} k) /dq x ^ for fc > u, we find that 

^ (£ n '°>°) =F,[gi,..., 9n _ z+ i], i = l,...,». 
Since £^°' 2 "+- = _ yPn+'O we obtain (putting j = n + 1 - i + a in (0EJ) 

^ ( £ n,0,2n+ CT > ) = _2^+„_ 4+1 ..,g ff+rl _ 4 ], i = l,...,n, 

(where as usual we denote q a — for a > n) where 

Vi™'" [<Zl, — , q a +n-i] — Fi[lU •") Qn — i+l] + /cr+ri-i+l (<Zl ! 9o-+n-i)- 

■ 

Due to their structure, equations 146|) make it possible to successively express (eliminate) the variables 
<7 CT+ i, . . . , q s+cr -i = q n as differential functions of qx, . . . , q a : 

q<7+l = fa+l [Ql>-- -)Qo] 

q*+2 = fa+2 9l)---)?<r , s 

. (49) 
qn = fn [qi, ■ ■ ■ , qa] ■ 

Let us first observe that performing the elimination (|49|l in the systems <(45ll must lead to cr-component 
systems of the form q tr — Z r [qi , . . . , q a ] , while for each system in 145(1 the last s — 1 components turn into 

some system of differential consequences of w^ 1 '^, . . . ^Wo™^ (and are zero on S Uo i.e. they are satisfied 
along any solutions of q t = Z r ). Therefore, after this elimination we obtain 

| q tr =X[qi,--- ,fo]> r = l,...,n (50) 

[ = if\. [wi , . . . , w a \ , i = a + 1, . . . , n 

with q = (g x , . . . , q a ) T and ip\ = q iitr - (Zy)\ 

Lemma 8 The first s vector fields Z°in l|5()(l commute: 

z i> z< j i>3 = fi- 



ll 



Proof. Obviously, in general, for — 1, . . . ,n, 



for some vector fields Vij that vanish on £ C M. only. Assume for a moment that for n 



cr—1 and 



for some i,j < s we have Vij 



(n,<r) 



(n,<r) 



/r, ir n - 7^ 0. As the vector fields Z i , were obtained by the 

reduction of the complete (in the sense of the infinite hierarchy) components of Z™, i?™,thus by increasing 
n — > n + [3, we do not change the form of Vy, which now has to be expressed by a higher dimension 



invariants u^ n+ ^' . Vij 



w 



(n+/3,cr) 



. ,W, 



(n+f3,a) 



But w 



J 3 



ir', '[qi,—, Qn+p] and lower 



dimensional invariants Wj 71 '^ are nonexpressible by the higher dimensional invariants w[ n+ ^' <7 \ so we get 
a contradiction. ■ 

We will now show that this procedure leads in fact to an infinite hierarchy of commuting flows. In 
order to do this, we will for a moment introduce a new index so that the vector fields in (|50|l will be 
denoted Z r ' as being obtained by reducing the n-component Killing systems (l4"Bl . 



Lemma 9 In the above notation 



Proof. According to l|42[l we have 



Z 



■n+l,o" 



for r = 1, . 



K+i+i =K+i tori = l,...,n-a. 

Thus, increasing s to s + 1 and keeping a unaltered (so that n changes to n + 1) the n — a equations (|49|l 
change to n — a + 1 equations 

q<y+i = f2+iki, -)9<r] = fa +i [qi, -,qa], i = i, ...,n- o- 

so that the variables q a +i, ■ ■ ■ , <7n are expressed by the same functions of q±, q a and a new elimination 
equation for q n +i appears. Moreover, (Z" +1 [<7])' 7 = (Z^lq]) 3 for j — I,..., a and r = l,...,s. Thus, 
replacing g CT+i by /™ +i [<ji, in (Z," +1 [gr]) J and in (^'[g])- 7 yields for j = 1, . . . , a and r = 1, .. . , s the 

same expression. But the first operation leads to the reduced vector field Z ' while the second - to 

X ,a - ■ 

Let us now take s + 1 instead of s (so that n — > n + 1) in l|45|l and 149(1 and perform the reduction. 
According to the above lemma we obtain the following sequence of s + 1 reduced systems: 



q tr 



■n+l, er 



ZJ for r 



1, . . . , s and q t 



Z 



■n+l,(T 
n+1 



i.e. we obtain the same sequence of s systems as before plus an additional system at the end of the 
sequence. Therefore, we see that this procedure leads to infinite hierarchies of commuting systems, since 
we can always increase n as much as we please without altering the already obtained systems generated 
in previous steps. 

The procedure described above can be generalized by using only some part of the equations in (|46|l in 
order to perform the elimination, since all of these equations are invariant along the flows of our Killing 
systems. Namely, we can skip the last a (with 0<a<n — a — 1 = s — 2) equations in l|46|) and use only 



the remaining equations (i.e. w a+l = I), ur_|_ 2 
Killing systems with n = s + cr + a — 1: 



0, 



IV 



(n,cr) 



0) to eliminate q a + a +i, ■ ■ ■ , 9n in the 



qu 



Z 



s+cr+Q — 1 



[qi 



■ q s 



i. 



l. 



(51) 



Thus, the index a indicates how many of the last equations in (|46|l we "forget" about. It turns out that 
the elimination that follows leads also to hierarchies of commuting equations. To see that, let us first 
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observe, that this elimination can formally be obtained by performing the above described procedure 
with the help of the Lagrangian density £ n .- a > 2n + cr + a ; since according to Theorem[S]we have 

Ei [C n - a ' 2n+a+a ) = E^ a (£™^ 2 «+-) = w ( ™ ] for i = a + 1, . . . , n. (52) 

Denoting Ei (jr^m^n+k^ as w ( n > m > k ) ^ w ] lere now ■ u ^™' 0,fc ) = w\ n (in the notation of ((46(1 1. the last 
tl — o — a Euler-Lagrange equations (invariants), associated with £™> _a ' 2 ™+< T + Q ; have therefore the form 



(n, — a,(7-\-a) (n,(r) (n.a) r i r, 

r a + a = w <y = w <y ' [Qi, -,9n] = 
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a+1 — ^ fa+ 



(53) 



Wn = w y n _ a ' = -2q a+a+1 + Lp y n _^[q u ...,q a+a \ = 0. 



These equations make it possible to successively express (eliminate) the variables q a + a +i, . . . , q n as dif- 
ferential functions of qi, . . . , q a + a , which yields 



Qa+a+l — <?cr+Q+l • • • j 1cr+a\ 
qn = qn [qi, ■ ■ ■ , qa+a] ■ 

Therefore, after this elimination the Killing equations ((511) take the form 



(54) 



1t r = M - 

j (n, — a,<T-\-a) (n, — a,(T-\-a) 

= <Pr W l >•••)<+« 



i = a + a + 1, . 



s + a + a - 1 (55) 



with q = (gi, . . . ,g CT+Q ) T and = qi. tr — (so that the reduced systems will have N = a + a 

components). Similarly as before, in Killing equations 1)51(1. only up to r = s the first er + a components 
are complete in the sense of the infinite Killing hierarchy 1(35(1 . As before, it stems from the fact that the 
first s vector fields Zl commute to zero: 



—a+a -=<J+a 



= 0, i,j = 1, . . . ,s. 



the proof of which is analogous as in the case a = but now we have to take n = s + a + a— 1. As 
previously, we can repeat the elimination procedure taking s + 1 instead of s (so that n increases to n + 1 
and k = 2n + a + a increases to 2(n + 1) + a + a = k + 2 while a and a are kept unaltered). By the 
same argument as before, this new procedure (with n + 1 instead of n) will lead to a sequence of s + 1 
autonomous (a + a)-component systems in which the first s systems will coincide with the corresponding 
systems obtained from the original procedure (with n). Thus, again we will obtain infinite hierarchies of 
soliton systems. 

6.2 Elimination for negative potentials 

We now present the second possibility of elimination - with the use of negative (rational) separable 
potentials. Again, our aim is to produce s (s € N) commuting a-component (a <E N) vector fields 
(evolutionary systems) from 1(34(1 and 1(21(1 . This time however we have to choose n = s + 2cr — 1 and the 
Lagrangian density j n order to create an infinite hierarchy of commuting flows. 



Lemma 10 The first n — a invariant equations \21\) with a > n ) i.e. with m = n — a and k = —n, 

have the form 

v[ n ^^-± +1 [ n -°\ qi ,...,q„}=0, (56) 

(n,a) 2<7„_i_|_i (n.a) r 1 n ■ o 

v\ ' = 3 h7i m, ■ ■ ■ ,q<?-i+i,q n -i+2, ■ ■ ■ >2nJ = 0, 1 = 2, ...,n-a. 

Qn 

where we denote, to shorten the notation, Ei (£» ! ™- cr .-' ! ) as t ,(™' <T ) anc [ q a — q w hen a < 1. 
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Proof. From the recursion l|27|l it follows that 

V}~ j) = Vt i] (Qn-j+u ...,q n ), j = !■■-, n (57) 
From this and from (|B5|I . (gJJ and (|B7|I . we have 

E i (L^-^) = G i [q 1 ,...,q„- i+1 ], i=l, 

E i (C n ' n -' , > ) = Gi[q n - i+1+ir ,...,q n ], i = a + 1, ...,n - a. (58) 
Moreover, by using Lemma we find 

~ = s h gi{q n -i+2, q n ), i = 2, n - ct 

and 

dq a+1 dq x ql ' 

Plugging all this into E t ^£, n ' n ~ a '~ n ) , i = 1, n - a, we obtain JSSJ where % [q] = G t [q] - g t (q). ■ 
Let us now consider the following Killing systems 

q, r = Z"[q], r = a + 1, . . . ,a + s, with n = s + 2a - 1 (59) 

We can use the n — a equations l|56|) to successively express (eliminate) the variables q a +\i ■ ■ ■ ,q n as 
differential functions of q\, . . . , q a . This leads to the elimination relations of the form 

q*+i = f"+i [9i> ■ ■ ■ A<A , i = l,...,n-a. (60) 
Performing the elimination l|60|) in l|59|l we obtain an autonomous sequence of s evolution equations 

qt r =K Wi>---,Q<t] , r = a + l,...,a + s (61) 

such that the vector fields Z a r mutually commute to zero. One proves this by the same arguments as in 
the positive case, since the first a components of all the vector fields in l|59(l are complete in the sense of 
infinite hierarchy (|35Jl ). 

Analogously to the positive case, we will now show that this procedure leads to an infinite hierarchy. 
As in the positive case, we will for a moment introduce a new index so that the vector fields in l|61|l will 
be denoted Z r ' as being obtained by reducing the n-component Killing systems (|59|) . 

Lemma 11 In the notation as above 

Z r+1 = Z r for r = g + 1, . . . , a + s 
Proof. Let us observe that, according to results lBT|l . ifHf and l|B9|l 

for i = 1, . . . ,n — a. (62) 

qj— ►qj + i, j=n-i+l,...,n 



V; — V) 



Thus, increasing s to s + 1 and keeping a unaltered (so that n changes ton+1) the n — a equations (|6U|I 
change to n — er + 1 equations 

Qa+i = .Ctl[qi,---,q<j} , qa+i+i = .C+Z+ifei, -Aa\ = f2+ilqi, -,q<A, i = l,-,n-a. (63) 

Observe that the last n— a equations l|6U|) express now (fc+i+i (instead of q a +i) as f™ +i [qi, qa]- On the 
other hand, due to (|34|) . by changing qi — > qi+i for all i > a in the sequence qt r = Z™ +1 [q] we transform it 
so that (Z; l+1 [q}Y -> (Z%l[q]) j for j = 1, . . . , a and r = a+ 1, . . . , a + s. Thus, inserting /™ +i [<?i, q a ] 
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instead of q a +i+i in (for j = 1, . . . , a and r = a + 1, . . . ,a + s) yields the same expression as 

inserting the same function f£ +i [qi, ■ q a ] instead of q a+i in (Z™ +1 [q]) J . But the first operation leads to 

the reduced vector field "Z^y" while the second - to ■ 

Let us now take s + 1 instead of s (so that n — > n + 1) in 1|59[1 and i|60|) and perform the reduction. 
According to Lemma ll II we obtain the following sequence of s + 1 reduced systems: 

9i CT + 1 = Z a+1 , q tr+1 = Z r+1 = Z r tOT V = <T + 1, . . . , (T + S 

i.e. we obtain the same sequence of s systems as before but shifted and an additional system in the 
beginning of the sequence. This first system can therefore be treated as a next system in some infinite, 
commuting hierarchy of vector fields. 

Let us also observe that we could use <|60|) to eliminate variables in Killing systems of the form (|45|l 
and this would lead to a system of s commuting evolutionary systems. However, this choice does not 
lead to any hierarchy: by increasing s to s + 1 we obtain a different sequence of systems. 

As before, this procedure can be generalized: we can use the first n— a— a equations (0 < a < n — a — l) 
in i|56|) to eliminate q a + a +i, ■ ■ ■ ,Qn from the following sequence of s Killing systems 

qt r =Z?[q], r = a + a + 1, . . . ,a + a + s with n = s + 2(a + a) - 1. (64) 

This elimination leads - similarly as above - to s commuting to zero N = a + a -component systems 
q t = Z a r +a [gi , . . . , q a +a] and by increasing s by 1 we always obtain a new system of the hierarchy at the 
beginning of the sequence. 

Next section contains some examples of the above described elimination procedures. 



7 Examples 

7.1 Elimination with positive potentials 

Below we will present some examples performed with the help of the (generalized) procedure described 
in the previous section. Soliton hierarchies are now classified by pairs (a, a), a = 1,2,..., a = 0,1,..., 
where N = a + a is a number of components in the systems of a given hierarchy. Assume we would 
like to construct first s members of the hierarchy. We have then to fix n = s + a + a — 1 and take 
first s Killing equations in l|51l) . Then, we have to eliminate coordinates q a +a+i, 9<t+q+s-i = 9n using 
invariants = 0, ...,w^f2 = 0- According to (|52*|l . these invariants can be generated, for example, 

from C n ' 0,2n+a by taking the equations E t (£™A 2 "+ <T ) — o, i = a + 1, ...,n — a. After the elimination 
procedure, soliton equations are represented by first N = a + a components of first s reduced Killing 
equations. Observe, that in this procedure the first soliton equation has always the trivial form q t = q x , 

q=(qi,---,q*+a) ■ 

Let us start with a one-field hierarchy: N = <j+a = 1. There is only one possibility here: a = 1, a = 0. 
We present how to produce first s = 3 flows which will be recognized as the first members of the KdV 
hierarchy. We have therefore to take n = 3 and k — 7. Killing systems (|34|l have the form: 

= Z\ 

q2,x 

q3,x + q\q2,x - q%qi, x ) = z \ ( 65 ) 
qiqs.x - qmi,x 





3 



q3,x 

qiq-i,x - 9392,x 



while the Lagrangian £ 3 '°> 7 is 



✓•3,0,7 I22 1 2 1 1 1 2 n , o 2 , o 2 a 3 1 5 

£ ' = 9m ~ 2 q2<ll > x ~ 11 r<l - r ' ' + 2l2,x - 2 9293 + 3^ q 3 + 3qiq 2 - 4q x q 2 + q x 
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and the Euler-Lagrange equations (|46|l for the above Lagrangian attain the form 

t<>2 = -2q 3 + 69192 - 4gf + \q\ x + qiqi, xx - q 2 ,xx = 
^3 3,1) = -2q 2 + 3qf - q x>xx = 0. 

These equations can be solved with respect to 92, 93 yielding (|49() of the form 

92 = -\qi,xx + \q\ , 93 = \qi,xxxx - §9191,** - \q\^ x + f 9i- (66) 

Substituting it to the above Killing systems gives (|3T7|> that read now explicitly as 

Si,ti = qi,x = Z\ { qi.t 2 = —\q\,xxx + 3qiqi,x = z 2 

= < 0=0 

= [ Q = \w XjX 

9x,t s — \q\,xxxxx — 5qi,xqi,xx — ^qiqi,xxx + ^q\q\.x — Z^ 
= \w\.x 

= -jwi.xxx + ^qiwi.x 

so that the first components qi tti — z\[qi\ are the first three flows of the KdV hierarchy while the 
remaining equations are just differential consequences of w±, which of course vanish on any S Uo . By 
taking larger s we can produce an arbitrary number of flows from the KdV hierarchy. 

Next, let us consider two-field systems: N = a + a = 2. There are two possibilities: (a, a) — (2, 0) 
and {a, a) — (1,1)- Therefore, as a second example we consider the case (a, a) = (2,0), and s = 3. 
We have now to take n = s + <j + a — 1 = 4 and k = 2n + a = 10. The Euler-Lagrange equations 
E 4 (£ 4 - ' 10 ) = w { t 2) = and E 3 {£ 4 - ' 10 ) = w { 3 ' 2) = can be solved with respect to q 3 ,q 4 yielding $® 
of the form 

93 = - j9i,xx + 3<ji<72 - 2q\ 

qi = \ql, x - \q2,xx - qm,xx - \q\ + 3qlq2 + \q\- 

Substituting it to two first components of Killing equations Z\ [q] , Z 3 [q] yields two first nontrivial members 
of another two-field soliton hierarchy: 

91, t 2 = q2,x /g^N 

92, (2 = -\qi,xxx + 4gig 2 ,x + 2q 2 q ltX - 6qfqi tX ' 

and 

91, t 3 = - \q\,xxx + 39i9 2 ,x + 39291, x - §q\q\, x 

92, t 3 = ~ 2 q2,xxx - f 9l9l,xxa; + 39 2 92,x + 691 9291,2; + 69? 92,1c - 189^91, x . 

In the second two-field case (er, a) = (1, 1), if we keep s = 3 unchanged, we have to take n = s+a+a—1 = 4 
and k = 2n + a = 9. From the Euler-Lagrange equations (|46|1 for C 4 ' 0,9 

E 3 (C 4 ' ' 9 ) = 4 4,1) = 0, E 2 (£ 4 >°' 9 ) = «4 4ll > = 
we can eliminate 93 and 94, which yields Q54JI ■ Explicitly, we obtain: 

93 = \q\,xxxx — \q\q\,xx — \q\, x + f 9i 

94 = \q2,xxxx — \q\q\,xxxx — \ q\, x qx, xxx — \q\, X x ~ q2q\,xx — § 91,2:92, x 

+ Mfqi^xx - \q\q2.xx + TqiQi, x + 3qfq 2 - \q\ + ^q 2 . 
Then, two first components of the Killing equations qt 2 — Z 2 [q] , q t3 = Z 3 [q] turn into 

91, t 2 — 92, z 

92, t 2 = -\q2,xxx + \q\qi,xxx + qi, x qi,xx + 49192,2 + 2q 2 9i,x - 6q 2 qi,x ' 
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and 

<7l,t 3 = — nQ2,xxx + oQlQl,xxx + Ql,xQl,xx + 39l92,x + 3(7251.2; ~ QQiILx 

— 1 _ 1 _ _ 7 _ 7 

92, £3 — 4 Q2.xxxxx ^QlQl,xxxxx Ql,xQl,xxxx ^Ql.xxQl,xxx Q2Ql,xxx 2 92,2?9l,2;2; 

+ ^QlQl,xxx — 3qiq2,xxx — \<l\,xl2,xx + ^-QlQl,xQ\,xx + 6gi q2l\,x + 3(j292,x 

+69192,2- + - 18gfgi, x . 

Finally, we shortly mention the three-field case: N = a + a = 3. There are three different hierarchies 
with the following first nontrivial member of each hierarchy: 
for (a, a) = (3, 0) : 

qi,t 2 = q2,x 

q2,t 2 = qz.x + qiq2,x - q2qi,x 

93, t 2 = -\q2,xxx - 12qiq 2 qi tX - Gqlq 2 .x + 39292,2- + 1<fy<l\& + 49193,2; + I0qfqi, x , 
for (a, a) = (2,1): 

qi,t 2 — q2.x 

qi,t 2 = qz.x + q\q2,x - QMi.a; 

93,*2 — —\q2,xxx + \q\q\,xxx + Q\,xQ\,xx — 129l929l,2- _ 6<7i<72, x + 3(7292,2; 

+ 29391,2; + ±q\q?,,x + 105? 

and for (cr + a) = (1, 2) : 

91, t 2 = 92,2- 

92, t 2 = 93,2; + 9i92,2- - 9291,2- 

93, t 2 — —\qz,xxx + \qiq2,xxx + \q2q\,xxx — \q\q\.xxx — lq\q\, x q\,xx + q\,xq%,xx + q2.xq1.xx 

- \q\,x + 29391,2; + ^qiq^.x - 12919291,2; - 69192,2; + 39292,2- + 109191,2- 

In general, for a fixed N = a + a, this procedure leads to TV different N— component hierarchies of soliton 
systems. As the field representation of constructed hierarchies is non-standard, it is not easy to recognize 
which hierarchies are known and which are new. We immediately recognized the KdV hierarchy. We also 
found that two-field hierarchy starting from (|(j7[) turns after the transformation 

til = — 3<7i + 292, u-2 = 2qi, x — > \ 2ix, t — ► \ 2~it 

into the 2-component coupled KdV hierarchy in the representation of Fordy and Antonowicz 23 . For 
example, the first flow of this hierarchy l|t)7|) turns into 



Ul.ti = jU2.xxx + \u 2 Ui, x + UiU 2 ,x 
U2M = Ui )a! + %U 2 U 2 ,x- 

(yielding (3.18) in [2Sj)- 

7.2 Elimination with negative potentials 

We start by presenting the first two (s = 2) flows of the only N = 1-component hierarchy that can be 
obtained within our scheme by using the negative separable potentials. Since N = 1 = a + cr, the only 
choice is to put a = 1, a = 0, which yields n — s + 2(a + a) — 1 = 3. The Euler-Lagrange equations (|56|l 
for 

2 2 

f*n,n— (t, — n /^3,2.— 3 *~ J2 ! " ' ' 

z - — *~ — o9l,2; 



•m,n-<r,-n _ ,.3,2,-3 _ * 2 _ 92,2-93,2- 9 3 ,2;92 _ 9l 92 

2" ' 93 293 9 3 2 9l 



attain the form: 

-1 - 93 9i,x2; = 0, 492 + 2qlq3,xx ~ 9393 2; = (68) 
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which allows for expressing q 2 and q 3 as differential functions of q\. 

93 = 93 [gi] = {-qi,x X y 1/2 

92 = 92 [gi] = -33 (53i l3aC s ~ l9:..r.r-9i., ,., . ) (-9i,x: C )~ 7/2 

(here and in what follows we only consider the positive solution for q n , otherwise we can change t — ^ —t). 
Substituting these expressions to the Killing systems (|64[) and performing the necessary derivations we 
obtain the following two commuting flows: 

9i,t 2 = (gafeiDs > 9i,t 3 = fefeiDa 

with the differential functions 92 [91] and 93 [91] given as above. After substitution u = —qx tXX the second 
equation turns into the well known Harry Dym equation while the first one becomes the second member 
of the hierarchy. If we want to produce a next member of this hierarchy we have to take s = 3. According 
to the general remarks in the previous section this new system will appear as the first system in our 
sequence of systems. 

Let us now consider two-field systems: N = a + a = 2. As before, we have now two choices: 
(a, a) = (1, 1) and (a, a) = (2, 0). We start with (a, a) = (2, 0). We have now n = s + 2cr — 1 = 5 and 
thus we consider the Lagrangian £p,n—cr,-n _ £5,3,-5^ -j^g assoc j a ^ ec j Euler-Lagrange equations H56fl can 
be written as 

95 (li,x + Zqi.xxQi - 1q 2 , xx ) -2 = 0, 2q 4 - qlqi. xx = 0, 
49395 - 95,^95 - 6 94 + 2 9s95,xx = 

and they can be solved to 

95 = /|[9i,92]=2 W - 1 / 2 

94 = /I [91, 92] = Mi :XX w~ 3/2 

93 = / 3 5 [9i, 92] = (-§«£ + ^ql, xx w + 2ww xx ) w~ 7/2 

where w = 2q\ x + 4:qi iXX qi — Aq2 >xx . Substituting it into (|64|l we arrive at the following two commuting 
two-component systems: 



91, t 3 = (/ 3 5 [9i,92]) K 

92, t 3 = 9i (fl[qi,q2]) x - (/ 3 5 [9i,92]) q\, x + (/I [9i, 92]) 3 



and 



qi,u = (f![quq2]) x (69) 

92,t 4 = 9i (/|[9i,92]) 2 . - (/f[9i,92]) qi, x + (/I [?i> ffa])^ 
The system (|69ll can be written more explicitly as 

qi M = 2 (2qi :XXX w - 3qi, xx w x ) w" 5/2 

92, u = (Miqi,x XX w - 6qiqi, xx w x - 4qi tX qi tXX w - ww x ) u>~ 5/2 . 

Finally, let us consider the case (<7, a) = (1,1). Again, we have n = 5, but this time we consider the 
Lagrangian £ n . n -< 7 - a *~ n + a = £ 5 > 3 >~ 4 . Its first n — a — a = 3 Euler-Lagrange equations 



= -qjqi, xx - 1 
= 4g4 + 2q%q 5tXX - (75175 

= 2(73(75 - 9594,x95,x + 949595,x - 3 94 + 9594, xx - 959495,xx 
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yield the following elimination equations: 



95 = fi [Ql\ = (-Qi,xx) 1/2 

94 = fi[qi] = -jz ( 5 9i,in - Mi,xxqi,xxxx) (-qi,x x y 7/2 (70) 

93 = fiiqi] = w p \9i] (-9i,xx)" 13/2 

where P[qi] is some complicated differential polynomial of q\ (homogeneous of degree 4 and of order 6) 
with integer coefficients. Substituting (|TU1) into the Killing systems we arrive at the following two 
commuting two-component flows: 



9i,U = (/I fed)* 

q2,t 3 = qi (/fML - fl[qi]qi,x + (/|[<Ji]) 3 



and 

p5 



9i,t4 = (/|[9l]) B 

qiu = qi (f![qi)) x - f![qi]qi,x + {f![qi}) 3 

The last vector field can be written more explicitly as 



qiM = W (~^0ww x w xx + 35wl + 8w 2 w xxx ) w 9/2 

92, u = (lQqi, x ww 2 x - 8qx <x w 2 w xx - 40qiww x w xx + 35qiw x + 8qiw 2 w xxx ) w~ 9/2 

where w = —qi tXX . Let us notice that in this case we obtain a hierarchy of systems such that every system 
is driven by its first equation which is a consecutive equation of Harry Dym hierarchy. One can see that, 
contrary to the positive case, if a > the obtained systems are always driven by its first a components 
that coincide with the corresponding systems from a — hierarchy. 



8 Conclusions 

In this paper we developed a method of unified constructing of Stackel systems and soliton hierarchies 
from the same common denominator in the form of separation relations (JJJ. We developed our theory 
starting from separation relations generated by separation curves of the form 

ffiA 01 +... + H n \^ = iA"> 2 + A fe , ft,n€N, m,keZ. (71) 

We performed a detailed, systematic construction of soliton hierarchies for the Benenti class of separation 
relations, given by the particular form of (|71|l . namely 

ffiA"- 1 + H 2 \ n - 2 +-.- + H n = ^A'V + X k . 

The results we obtained are hopefully only a first step of a new research program. The next step of 
this program would be finding out a way for systematic constructing of other soliton hierarchies from 
different classes of separation curves l|71ll . when the sequence (/3i, ...,f3 n ) differs from (n — 1, ...,0). The 
next - nontrivial - step would be to extend the theory to the case of polynomial separation curves © with 
(«!, ...,a n ) 7^ (0, ...,0). We expect by presented procedure to generate not only the majority of known 
soliton systems but also to construct in a systematic way a vast number of new integrable hierarchies. 
One should also investigate the possibility of " prolongation" of standard integrable structures of separable 
systems (such as integrals of motion, bi-Hamiltonian structure) onto the corresponding evolutionary 
hierarchies of PDE's. 
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9 Appendix A 

The involutivity of jj( m,fe ) and Hi m ' k ^ (I15fl leads to the following relations imposed on (A), (A), (A) 
and T/ r (fc) (A) [Hj: 

0, z = l,...,n, (Al) 



dk 

dv 

5A,; 



= TTT^-T. ^ fc ' all /7i, r (A2) 



—!—=vl— i— for alii, r and A. (A3) 
CMi aAi 

We will prove here the relation pOjl . i.e. 

ELx^ (£("< m ^) X i>tr - Eti^ (4 n ' m ' fe) ) A*,. - 0. (A4) 

First, let us consider the geodesic case. Due to ljT5t we have Ai,t r = v l r \^ x where are eigenvalues of 
if r (see H161 ). For the geodesic Hamiltonians 

£ (n, m ,0) = l^ 5 M A 2 K and £(^0) = l^ 9 ( r ) v i A ^ 

i=l i=l 

so that in this case the left hand side of l|A4[) attains the form 

££=i#< (£ (n ' m ' 0) ) A i>tr - EHi^i (4" ,m,0) ) Ai >a 



iE|-fc.) I -;v-E5( ! r i >,).iv 

i,k=l i—1 

— 1 2,/C — 1 

n , n j z \ 



dx V" ""J r 4 < K cte 

1 i=l 



1 n ^ (m) I n a fc n 

1 "5fcfc /\ \2 / j k\ \ I 1 (m)C w r /> >2, (m),2 



dx 



since expression in the last parenthesis equals to 

2 5 fefe ~~ 2^ 5 « ^QaT ' x 

i.k— 1 1 



2 ^ <9A* v ' ' ' 2 9A 

i,/c— 1 1 



«\. ~ V W r w rJ 



where the last equality follows from formula (|A2|I which can be written in equivalent form 

d\i dk 

(that is in fact also valid for k = i). Thus, the statement has been proved for geodesic densities. For the 
potential parts: 

'dv} k] l dV r (ky 



ruE, (vf >) a.,«, - (*•) a,.. - g ( - ^- ) A,, - 

due to (|A3|) . This concludes the proof. 
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10 Appendix B 

We will prove here Theorem^! The relation (|40(l is a consequence of ()37|) and Ij39(l of LemmaEl Indeed, 
by 

a(£ „,„ t) _ 1 £ <5^__ (9()i (j; • („).) . 



and 



jj | y(2n-W-g— i-i+ct) j lcmma|51 El ^£rum,k^ 



dx , 

\i=l / 

The relation l|41|l is just a rewritten form of (??). 

Since in what follows we will compare separable potentials with different n, in the rest of the proof 
we will use temporary extended notation for potentials in the form Vi"' . From (|28[1 and (|47|l it follows 
that 

v n,(n+k) = yn+Un+l+k)^ fe = _^ ^ _ ^ ^ 

We prove now (|42H . Using the relation (|37|) for r = 1 we obtain 

i n m/ n '( 2 ™-*-i) A ( n \ 

t,j=i Vi— l / 

- 5 S ^ («). («)• - £ IX*-" (•).) (B2) 

— l \i— 1 / 

and in a similar way we have 

i,j'=l \i=l / 

= 2 E ^ (*). (»). - Tx 1 E^ (*)- 

ij = l \i=l / 

tnn i ^ a^n-w-i+i) d / » \ 
= 2E ^ (»>.(*). -5; (E^i '(»). • 

i,j=l 1 \i=l / 

The equality (*) is due to he fact that ]/™+ 1 '( 2 " _4_ '+ 1 ) = for i = n+l and similarly y^+ 1 '(- 2n - t -i— l + 2 ) 
does not depend on qi for i = n + 1 or j = n + 1 (this follows from (|42[) and H57|) ~). Thus 

E t (£ n <°>°) = E l+l (£« +1 - '°) , Z = l,...,n. (B3) 

Moreover, from (|3(J[1 it follows that 



V-n,(2n+cr) T/ n,(2n+(T-l) T rn,(n+<r) 

1 — - ... - (fol^ , 



hence, for I > a 



oqi dqi dqi 



dv n,(2n+*-l) g n,(n +g - i+ l) n+ ^_ 
"?1 5 --Q-n 5 ^1 
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On the other hand we have 

T/ n+l.(2ri+o-+2) _ T/ n+l,(2n+<7+l) T/ n+l,(n+CT+2) T/ n+l,(n+CT+l) 

V\ — —QlVi — ■■■ — QnV 1 — q n +l v x , 

hence, for I > a and according to IjBlf) and l|47|) 

o T rn+l,(2n+o-+2) QTr n+l,(2n+o-+l) o Tr n+l,(n+o-+2) 

av l av l av l T ,n+l,(2n+er-i+l) 

o = 5 ••• - 3" ^ 

Oqi+i dqi+i 

njrn+l.(2n+a+l-l) aT ,ri + l.(?i+<r+2-i) 

- 2ll „ dV l T/ n+l,(2n+<7-;+l) 

= -qi ~ ••■ - q n 3 V-y 

dqi dqi 

_ „ dV l „ dV l T/ n,(2n+<7-0 

dqi dqi 



dv; 



n,(2n+er) 



dqi 

and from (|37fl it follows that 



gyn,(n+s) Qyn+l,(n+s+2) 

— ^ = — ^ , 0<s-Z + l<n. 



-11 +1 



So, from l|B3|l and l|B4|l for a < I < n equation lj4"2"|) is fulfilled. 
Now, we pass to the proof of relations l|43|) . First, we have 



1 » o T/ n,(n+o— i-j-i+l) A I n , 



On the other hand 



Ei = £ £ _ {qi)x {qj)x -^(Y. Vr Un+ — l+1) (ft). 



By (|42|) and l|57|) . for Z < cr the last term in both sums does not contribute. Moreover, according to 
and the fact that y"^ 2 ") — q n +i = V^ +1 '^ 2n+1 \ we have 

n,(n+fc) ot ^n+l,(n+fc+l) 



) " ' - J ' **J 



9gi 9gi 
hence 

1 n Q T /n,(n+o— i-J-i+l) A ( n 

Ei (£ n +1 ,„ +1 -., 0) = 1 £ ^ ^(n+.r-i-O (ft)a 

= (£"."-".°) . 

Finally, we prove the relation (|44|) . From the negative recursion i|27|) . we have 



VT fc) (g„_fc + i,...,g„) = K; (g„_ fe+ 2, ■ ■ ■ >q n +i), k = l,...,n. 

Qi — >Qi+l, i—n—k+l,...,n 



r n+X,{-k) , 
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Then 



^ - \ E ^^Ss^ o. (*). - 1 E H*^ o.) 

i,j=i T \t=i / 

1 ™ m/™>(2"-»-:M+l) ri ( n , x \ 

- 5 E ^ — (*>, <•>. - s E vr'-'-" fe), (bt) 

— 1 \z— 1 / 



and 



K+i {£n +l, n +l-„ 0) = 1 ^ ( ^ ( ^ _ d_ ^ ^.(n^+D 

— 1 y +0" \i=l / 

5 E o. w. - s E («>.) 

i,j'=l \i=l / 

1 " gyn+i,(2 n -i-j-i+i) d ( n N 

2 E o 1 dqn+1 fe+o, - s (E fe + o^ 



As for i,j = there is no contribution to the sum, so according to l)B6)l we have 

E a+l {C n > n -^\ qo ^ qo+1 = E a+l W) , l = l,..., n -a, j = l,...,n. (B8) 

From $57$ and (|B6|> we have 



dv; 1 



^ - 1 -,i=l,...,n. (B9) 



%+l dq. 

9j— > 9i+l, 3=1, ...,-n 



that together with l)B8)l proves the relation 1)31)1. 
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